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We discuss progress in understanding the light and heavy quark excited hadron 
spectrum from Coulomb gauge QCD. For light quark systems we highlight the 
insensitivity to spontaneous chiral symmetry breaking, which predicts Wigner 
parity-degeneracy in the highly excited hadron spectrum and allows the quark 
mass momentum dependence to be experimentally probed. For heavy quark 
meson decays we invoke the Franck-Condon principle, a consequence of small 
velocity changes for heavy quarks, to extract qualitative, but model indepen- 
dent, structure insight from the momentum distribution of the decay products. 
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1. Motivation for Coulomb gauge QCD studies 

Recent years have witnessed steady progress in both understanding and ap- 
plying Quantum Chromodynamics in Coulomb gauge. Much of the attention 
has been formal, focusing on the such issues as the behavior of the Green's 
function in the deep infrared |T], the Slavnov- Taylor identities j2] necessary for 
a complete renormalization treatment in this gauge and the structure of the 
ground state vacuum [3]. In this paper the thrust is more phenomenological 
as we detail how an approximate Coulomb gauge formulation |3] can provide 
important hadron structure insight. While the QCD treatment is not rigorous, 
we have been able to identify several fundamental issues which a more precise 
Coulomb gauge QCD analysis can resolve. 

In addition to the advantages of Coulomb gauge QCD, it is worthwhile to 
also comment on this formulation's major limitation involving calculating ob- 
servables that connect different reference frames. Since the Coulomb condition, 
V • A" = 0, is specified for one chosen frame (typically the hadron center of 
momentum), the computation of a form factor or an electromagnetic transition 
matrix element, which requires hadron wavefunctions in two different frames, 
entails employing boost operators [5]. This introduces a formidable difficulty 
since in this gauge these operators are as complicated as the exact Hamiltonian. 
Although one can adapt methods developed for the relativistic quark model [5] , 
unambiguous, reliable results become very difficult to obtain. Related, fragmen- 
tation functions and parton distribution functions, more conveniently calculated 
in light front quantization, are also difhcult to formulate in the Coulomb gauge. 



In contrast, Coulomb gauge QCD is much more ideal for spectroscopy and 
hadron structure including the vacuum which is difficult for light cone formula- 
tions. Solving the equation of motions variationally either on the lattice or by 
approximate model truncation allows study of the entire spectrum with masses 
of arbitrary excitations (in the rest frame), any spin and all other quantum 
numbers. The variational treatment is possible because this gauge is physical 
(contains no spurious degrees of freedom or ghosts), involving only transverse 
gluons. Gauss's law, which is essential for confinement, is explicitly satisfied 
and incorporated as a constraint. However, the resulting interaction kernel is 
complicated. Fortunately, with insight provided by lattice QCD, it is amend- 
able to modeling using a more calculable confining kernel containing a linear 
potential as used in Ref. [3] and further supported by other approaches [?]• 
This approximate framework provides several new insights into the structure of 
highly excited mesons and baryons. One is the role of chiral symmetry for the 
complete hadron spectrum. As detailed in section 2, because chiral symmetry 
is incorporated (unlike constituent quark model formulations), one can clearly 
understand why the light quark, low mass pseudoscalar mesons are strongly 
governed by spontaneous chiral symmetry breaking while the high lying spec- 
trum is insensitive to chiral symmetry breaking leading to parity doublets. This 
is difficult to explain with other methods which generally do not respect chiral 
symmetry and/or can not predict the excited spectrum. For example, for heavy 
quarks non-relativistic QCD (NRQCD) [S] reasonably describes the low lying 
spectrum. However, it becomes inadequate for strongly excited states due to 
the attending large quark momenta. Again, the Coulomb gauge formulation 
naturally extends to systems involving arbitrary large momenta. Further, in 
addition to mass spectra, one can also learn about the momentum structure of 
the hadrons in strong decays, due to the large scale provided by the heavy quark 
mass and high momentum, rngv 3> Aqcd, which is discussed in section 3. This 
provides decay signatures which may be useful for identifying exotic hadrons 
and is analogous to the Franck-Condon principle of molecular physics. Section 
4 discusses new insights through future measurements at Jefferson Lab with the 
anticipated 12 GeV upgrade. Finally a summary of key results provided in this 
paper is presented in section 5. 



2. Insight into chiral symmetry 

Chiral symmetry is said to be hidden (or spontaneously broken) by the 
condensation of quarks in the ground state (vacuum) of Quantum Chromody- 
namics. This feature, documented in both Landau and Coulomb gauge studies, 
is expected to be of lesser importance in the excited spectrum as illustrated by 
a simple analogy. Consider two identical solid bridge supporting pillars in two 
different shallow ponds, one salt and the other fresh water (thus breaking the 
symmetry between them). Their net weight is essentially the same (approxi- 
mately symmetrical) even though one experiences a slightly larger buoyant force 
from the salt water. Two highly excited, but opposite parity, hadrons are anal- 
ogous to the two pillars since their quarks have (on average) large momenta but 
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now the running quark mass is small [5], so chiral symmetry is approximately 
restored. The first statement, that the momentum is larger for excited states, 
is supported by the relativistic virial theorem [10) and also explicit variational 
model computations for momentum distributions as detailed in our Coulomb 
gauge approach [TTl[Tlj. The second statement, that the running quark mass is 
small for large momentum, is a consequence of asymptotic freedom, well-known 
from early QCD studies |13) . 

An immediate consequence of chiral symmetry restoration is the reappear- 
ance of Wigner symmetry in the highly excited spectrum, with doublet partner 
states of equal mass but opposite parity [13]. For three quark configurations, 
relativistic chiral approaches like ours predict that states come in quartets, re- 
ducible representations that split into two doublets [15]. As an example, take 
the family of spin 1/2 excitations of the isospin 3/2 A. We have used the Cor- 
nell approach to Coulomb gauge QCD with harmonic oscillator interactions and 
a large, harmonic oscillator basis model space to diagonalize the Hamiltonian. 
All orthogonalizations and symmetrizations are performed numerically and the 
resulting A spectrum is plotted in Fig. [l] for the first 20 eigenvalues. Note that, 
while there are variations for the lower states, parity doubling is indeed a fea- 
ture of the higher spectrum illustrating chiral symmetry in the Wigner mode. 
Lattice gauge simulations and also Dyson-Schwinger approaches have difficul- 
ties generating excited state spectra which is easily computed in Coulomb gauge 
based models. 



Paiity doubling in the A spectrum? 
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Figure 1: Excited / = 3/2, J = 1/2 A spectrum. Parity doublets emerge, 
especially for higher excited states. 
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A potential problem investigating parity doublets experimentally, especially 
if the spectrum is numerous, is to identify the correct corresponding Wigner 
partners. This is where theoretical calculations, which do not have this ambi- 
guity, can be of great benefit as we again document in Fig. [2] where the doublet 
separation energy (splitting) is plotted for each J. The parity degeneracy clearly 
emerges for increasing excitation energy and J. This variational Monte Carlo 
investigation is still in progress and the attending error bars have not been 
shown. 
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Figure 2: Predicted baryon spectrum from variational Monte Carlo calculations 
(boxes) compared to data (solid circles) from Ref. [15] . 



Related, but perhaps more interesting, one can use the onset of this degener- 
acy to probe the quark mass momentum dependence (running) in the infrared. 
The quark mass in our approach appears both in the quark spinors and in the 
QCD Hamiltonian. For small quark mass and high momentum, a series expan- 
sion of the spinors and the Hamiltonian in the parameter m{k)/k reveals several 
insights. Expanding the spinors 
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with E{k) = a/Zc^ + m{k)'^. The first term is chirally invariant and the second 
chiral-symmetry breaking. Note that the lower component in each term has the 
opposite sign. Instead of directly expanding the QCD Hamiltonian, expand the 
Hamiltonian matrix elements computed in the Hilbert space spanning highly 
excited resonances, labeled n, where the average momentum, < fc >, is large 



i\HQ^'^\n) = {n\H^^'^\ 



i{k) 



(3) 



Again the first Hamiltonian term, labeled by xsi is chirally symmetric while the 
second, labeled xb, is not and involves nonchiral, spin-dependent potentials in 
the quark-quark interaction. 

Now recall for zero mass current quarks, the chiral charge, 



dx (x)75V'(x) 



(4) 



commutes with the QCD Hamiltonian. Nevertheless, chiral symmetry is still 
spontaneously broken by the ground state vacuum since QsjO) 7^ 0, leading to 
a large constituent quark mass in the quark propagator, pseudoscalar Gold- 
stone bosons and no parity-degeneracy in ground-state baryons. Insight about 
the chiral charge operator follows by performing a BCS transformation to a 
quasiparticle basis with operators B and D defined by 
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Then using the standard Fock operator normal mode representation for the 
quark field 



*W=E/(^["^W^Ac(k)+«A(-k)4c(-k)]e^''-6c , 



(7) 



and substituting for the bare spinors, the chiral charge can be expressed as 
dk ^ k 
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(<T . k)AA' (i3lc(k)SA'c(k) + i?l,c(-k)DAc(-k)) 
^ i^a2)^x' (k)^lc (-k) + Bye (k)i^AC (-k) 



The first two terms between the square brackets are quark and anti-quark num- 
ber operators, i.e. operators conserving the number of particles, flipping spin 
and parity. For m{k) << fc, it dominates the third and fourth terms represent- 
ing pion creation and annihilation. As argued in Ref. |12j . repeated action of 
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the chiral charge on a three-quark state produces a quartet of states, two of each 
parity, that dynamically split into two doublets of parity partners. Moreover, 
the mass splitting between partners is a direct measure of m(k) and they are 
degenerate when m{k) vanishes. 

In order to link this parity doublet mass splitting, AM = \M^^^ — M^^^ |, 
to the running quark mass, we examine the lowest-lying parity doublets for 
increasing spin J and incorporate the following key elements. 

• Regge trajectories, J = ao + aM^^ ''-^^ aM^'^ . 

• Relativistic virial theorem, (fc) — ^ C2M* — ^ ^VJ. 

• Canceling of the chirally invariant term, < n\H^g^\n >, in AM so that 
AM « Af± and AM < ^H^^'' > ^ 03=^^ < H^j^'' >. 

• In H'^g^, the spin-orbit term, L • S, is crucial to correct the angular 
momentum in the centrifugal barrier from to the chirally invariant 

+ 2L.S = — |. Due to the sign difference in the helicity-dependent 
term, —a ■ k, in the spinor, the spin-orbit term in H^g^ adds to the 
mass splitting instead of cancelling out as it does for H^g^. Since the 
centrifugal barrier scales like for high J, the spin-orbit term scales 
with one power of J less 

Combining these elements yields 

AM ^ ^^m{<k>)J-K (9) 

This equation linearly relates the parity splitting (an observable) to the running 
of the quark mass at the average momentum < k > for that splitting. Note 
that < k > will vary with splitting, increasing for higher excited states (and 
J). Although there are still unknown constants that are gauge dependent, a 
useful power-law running of m(< k >) (or equivalently m{k)) can be obtained 
model independently from data at high spin J by fitting the splitting to the form 
AM (X J^* and determining the exponent i. Then, the power-law behaviour of 
the running quark mass is given by 

m{k) cx fc-2*+2 . (10) 

As indicated in Fig. [2] the known masses of lowest-lying A resonances for J = 
1/2, ••• ,15/2 is currently not sufficient to derive this exponent. Establishing the 
masses of the parity doublets for spins J > 9/2 would resolve this and provide a 
running quark mass momentum dependence as schematically indicated in Fig. [3j 
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Figure 3: Cartoon showing the running quark mass as a function of momentum. 
The solid line is a power law whose exponent can be obtained from the decreasing 
parity doublet splittings for increasing J, as illustrated in Fig. [2j 

3. Franck- Condon principle and heavy hadron decays 

In general, calculating most hadron observables such as masses, widths, or 
partial branching ratios entails integrating over quark momenta which precludes 
direct, experimental study of momentum distributions. Here, we submit that it 
is possible to access this distribution by measuring the decay hadron momenta 
of open flavor, light-heavy quark qQ and qQ systems from highly excited heavy 
quarkonium QQ states. This assertion is based upon two points. The first is 
that QQ states with high excitation energies also have high momenta. The 
second is the heavy quark "velocity superselection rule" which states that in 
heavy quark meson decays the velocity of the heavy quark is conserved and not 
affected by the non-heavy degrees of freedom (light quarks and gluons). This 
selection rule is a common ingredient in heavy quark effective theories such as 
HQET and NRQCD and is due to the clear scale separation between the large 
mass of the heavy quark and Aqqd (as well as the running masses of light 
quarks). In particular, the NRQCD Lagrangian at leading order in ^qcd /^nq 
for heavy quarkonium is 

--(^•'^")("'^o-™^;^£,)(|:)-- 
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where the velocities, v and w, are the same for the quark creation and destruc- 
tion operators, /i^^^ . It is clear that the operator structure of this Lagrangian 
preserves the heavy quark velocity. Coulomb gauge QCD also takes this form 
in the heavy quark limit as can be seen from the explicit Hamiltonian |3] . 

Returning to accessing quark momentum distributions and invoking these 
two principles for a highly excited QQ decay to open flavor mesons, — >■ qQqQ, 
in leading order, the velocity (and momentum) of each heavy quark is the same 
before and after the decay. Then the quarkonium momentum distribution is 
approximately given by the momentum of the final state mesons containing 
the heavy quarks, which is directly measurable. This only applies to 3 or 4- 
body final states since 2-body final states have fixed momentum due to over 
all momentum conservation. Technically, in the limit of infinite heavy quark 
mass these principles, combined with momentum conservation, would forbid a 
two-body final state decay. This implies that for finite mass, decays to two- 
body states should be suppressed relative to three or more hadron final states. 
The above decay features should therefore be more pronounced in bottomonium 
transitions compared to charmonium decays. 

This effect [11] is very similar to the well-known Franck-Condon principle 
in molecular physics. It states that during an electronic transition, the nuclear 
degrees freedom (position and momentum) are essentially frozen due to the 
large scale difference between the electron and nuclei masses. Therefore, the 
wavefuncion associated with the motion of the nuclei remains the same after 
the transition. 

It is thus possible to reconstruct the excited state quarkononium wavefunc- 
tion by measuring the final state heavy meson momentum distribution. In 
particular, a highly excited pure QQ state has a definite structure of Sturm- 
Liouville zeroes (nodes) in the wave function, whereas a QQg hybrid state (that 
is, quarkonium plus a gluonic excitation) with the same mass, say in the ground 
state, would have no Sturm-Liouville zeroes (no nodes). Experiment can now 
distinguish between these two hadron structures. 

As an illustrative example consider the decay of T( 10860) to Bg Bg and 
B*g Bg final states. This parent state is generally accepted as an excited two- 
body quarkonium configuration, i.e. having several nodes. The momentum 
distribution for this state was calculated [TB] using a Coulomb gauge model 55 
wavefunction having four nodes and is plotted in Fig. [4] Also plotted (dotted 
line) is the same distribution but now a decay mechanism is included (^Pq scalar 
model) . Both distributions have also been multiplied by phase space factors and 
exhibit a clear node structure (depicting additional nodes requires a non-linear 
graph). Also plotted is the recent Belle Collaboration data [17] for T(10860) 
decay to Bg Bg and B*g Bg including 2, 3 and 4-body final meson states. Notice 
that there is reasonable agreement, especially the correspondence for the first 
Sturm-Liouville zero, affirming the validity of the Frank-Condon principle for 
heavy quark systems. 

Finally, one can also construct a different momentum signature to distiguish 
between QQg hybrid candidates and possible tetraquark states with similar 
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Figure 4: Decay widths for T(10860). Solid line: predicted momentum dis- 
tribution using the Frank-Condon principle for a 5S quarkonium wavefunction 
multiplied by phase space factors. Dotted line is the same as solid but now a 
decay mechanism is included. The dip in the data (dots) JJT^ is correlated with 
the first node in the momentum distribution. 



mass (both without Sturm-Liouville zeroes). Now it is necessary to study four- 
meson decays such as Z — ^ DDDD, or DKKD for which one can obtain, from 
the momentum distribution in the final state, the off-plane correlator, 

rr/ ^ l(PiXP2)-P3p 

n(pi,P2,P3,P4) = ^ , (12) 

n [ip»xp,i]V2 

1.3 = 1 

that reflects a distinctive non-planar structure from momenta combinations not 
present in 2-body or 3-body states. Consult Ref. [16] for a more complete dis- 
cussion on using this signature to identify tetraquark states. 

4- Future excited spectrum studies at Jefferson Lab 

As discussed at this conference, one of the goals for the 12 GeV upgrade 
of CEBAF is to explore the highly excited meson and baryon spectrum with 
existing experimental halls also having upgraded instrumentation. This will 
facilitate electroproduction measurements of excited baryon resonances with 
higher spins, J, allowing completion of the analysis involving Eq. |10[ Further, 
finding these excited states will bring closure to quark models which have long 
predicted their existence. There is also a QCD motivation, to understand the 
onset of chiral symmetry breaking in the low spectrum, by clearly linking m(/c) 
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from the asymptotically free regime to the eonstituent quark mass. This is 
possible because m/k is a small parameter in the excited spectrum. 

Equally, if not more, important is the construction of the new Hall D to 
conduct Glue-X experiments to search for exotic hadrons with ghionic degrees 
of freedom. This hall will also permit exploring the excited meson spectrum, 
complimenting the other halls in this endeavor. Since discovering exotica is a 
major goal of Gluc-X, having a constraint on momentum distributions would be 
very useful and this may be possible again using the Frank-Condon principle. 
Consider for example the (j){2nQ) observed in electron-positron collisions by 
BaBar in the reaction 

e-e+ (?i(2170) ^!>/o (980)7 • 

Jefferson Lab will be able to produce it via 7p ^(2170)p. Since this reso- 
nance has many open decay channels [K'^K^ , K^tttt, etc) measuring the 
momentum of the two charged kaons may be useful to gain insight regarding the 
momentum distributions of the strange quarks in (^(2170). This would permit 
testing the presumed 3S quark model assignment (two Sturm-Liouvillc zeroes), 
or even documenting possible exotic components. Although the strange quark 
mass, rus, is now comparable to the strong interaction strength, Aqcd, the 
scperation scale might still be sufficient for qualitative insight since the strange 
quark momentum, p = nisv/ ^1 — v"^ /c^, is large according to the relativistic 
virial theorem. 

5. Summary 

As discussed in the previous pages, the Coiilomb gaiigc formiilation of QCD 
has several advantages and provides an opportunity for dynamical insight re- 
garding hadron structure. While not readily amendable for form factor pre- 
dictions, it provides a natural framework for spectroscopic studies for highly 
excited states as well as exotic configurations involving explicit gluonic degrees 
of freedom. The role of chiral symmetry has been clarified, both for light quarks, 
where spontaneous symmetry breaking by condensates in the vacuum leads to 
a small mass pion, as well as highly excited states, which are insensitive to 
this symmetry, yielding parity doublets. Further, measuring the mass splitting 
between doublets allows determining the quark mass momentum dependence 
which connects the asymptotically free current mass to the much larger con- 
stituent quark value. Finally, it has been demonstrated how the Frank-Condon 
principle can provide information about the quark momentum distribution in 
heavy quarkonium from decay measurements, some to be performed with the 
Jefferson Lab 12 GeV upgrade. 
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